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ABSTRACT 


A  high-resolution,  multi-level,  primitive  equation  model  is  initialized  with 
climatological  data  to  study  the  combined  effects  of  wind  and  thermohaline  forcing  on  the 
ocean  circulation  of  the  California  Current  System  (CCS).  The  ocean  circulation  is 
generated  by  the  model  using  a  combination  of  climatological  wind  stress  and 
thermohaline  forcing.  In  the  first  experiment,  the  effects  of  thermohaline  forcing  alone  are 
evaluated,  in  the  second  experiment,  previously  conducted,  the  effects  of  wind  forcing  are 
isolated,  while  in  the  third  experiment,  the  combined  effects  of  wind  and  thermohaline 
forcing  are  looked  at.  The  results  from  the  combined  experiment  show  that  even  though 
the  effects  of  wind  forcing  dominate  the  CCS,  the  additional  effects  of  the  thermohaline 
forcing  results  in  the  following:  the  seasonal  development  of  a  poleward  surface  current 
and  an  equatorward  undercurrent  in  the  poleward  end  of  the  model  region;  an  onshore 
geostrophic  component,  which  results  in  a  temperature  front  and  stronger  surface  and 
subsurface  currents  between  Cape  Mendocino  and  Point  Arena;  and  a  region  of  maximum 
eddy  kinetic  energy  inshore  of  ~125°W  between  Cape  Mendocino  and  Point  Arena, 
associated  with  the  temperature  front.  These  model  simulations  are  qualitatively  similar 
to  recent  hydrographic,  altimetric,  drifter,  and  moored  observations  of  the  CCS. 
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I.  INTRODUCTION 


Three  types  of  water  masses  are  found  in  the  California  Current  System  (CCS): 
Pacific  sub- Arctic,  North  Pacific  Central,  and  Southern  waters  (Hickey,  1997).  The 
Pacific  sub-Arctic  water,  characterized  by  low  salinity  and  temperature  and  high  oxygen 
and  nutrients,  is  advected  equatorward  in  the  CCS.  North  Pacific  Central  water, 
characterized  by  high  salinity  and  temperature  and  low  oxygen  and  nutrients,  enters  the 
CCS  from  the  west.  Southern  water  characterized  by  high  salinity,  temperature  and 
nutrients,  and  low  oxygen,  enters  the  CCS  from  the  south  with  the  poleward  undercurrent 
(Lynn  and  Simpson,  1987). 

The  climatological  mean  CCS  is  a  classical  eastern  boundary  current  (EEC) 
system,  which  consists  of  several  large-scale  currents  (see  Figure  1).  The  predominant 
flow  is  the  California  Current  (CC),  which,  in  the  mean,  is  a  broad  (-1000  km),  relatively 
slow  (-10-30  cm/s),  equatorward  surface  flow.  It  flows  year-round,  and  extends  to  -500 
m  depth.  The  second  basic  flow  is  the  California  Undercurrent  (CUC),  which  is  a 
narrower  (-10-40  km),  relatively  weak  (-2-10  cm/s),  poleward  subsurface  flow.  It  can 
vaiy  seasonally  and  is  strongest  at  -300  m  depth.  The  third  flow  component  is  the  Inshore 
Countercurrent  (IC),  which  is  known  as  the  Davidson  Current  (DC)  north  of  Point 
Conception.  This  is  also  a  relatively  weak  (-5  cm/s)  poleward  flow,  found  during  the  fall 
and  winter  at  the  surface  and  near  the  coast.  There  are  other  ICs  known  as  the  Southern 
California  Countercurrent  (SCC)  to  the  south,  and  the  Southern  California  Eddy  (SCE) 
inshore  of  the  Channel  Islands  within  the  California  Bight  (Hickey,  1979, 1997). 

Recent  observational  studies  have  shown  that  the  CCS  is  not  the  quiescent,  stable 
system  of  currents  with  a  well-defined,  unchanging  structure  suggested  by  Figure  1. 
Rather,  the  flow  fluctuates  greatly  in  both  time  and  space  (Chelton,  1984).  There  exist 
mesoscale  meanders,  eddies,  filaments  and  jet-like  surface  currents,  which  are 
superimposed  on  the  large-scale  flow  (e.g.  Bernstein  et  al,  1977;  Stmb  et  al,  1991).  The 
combination  of  these  features  has  led  to  a  new  conceptualization  of  the  CCS  as  a  system 
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of  currents  with  filamented  jets  and  mesoscale  eddies  modifying  the  mean  flow  (Mooers 
and  Robinson,  1984). 

Recent  numerical  modeling  results  from  Batteen  (1997)  have  shown  that  both 
wind  forcing  and  coastline  irregularities  are  important  mechanisms  for  the  generation  of 
many  of  the  observed  features  of  the  CCS.  In  particular,  the  alongshore  component  of  the 
wind  stress  has  been  shown  to  be  a  key  ingredient  for  generating  realistic  vertical  and 
horizontal  structures  for  the  surface  equatorward  and  subsurface  poleward  currents.  With 
such  structures  the  currents  are  baroclinically  and  barotropically  unstable,  resulting  in  the 
generation  of  meanders,  filaments,  and  eddies.  Irregularities  in  the  coastline  geometry 
have  also  been  shown  to  be  important  for  “anchoring”  upwelling  and  filaments  as  well  as 
enhancing  the  growth  of  meanders  and  eddy. 

Recent  dynamic  height  analyses  by  Batteen  et  al.  (1995)  have  shown  that  the 
distribution  of  both  temperature  and  salinity  in  the  CCS  can  be  important  in  defining  the 
large-scale  circulation.  In  particular,  while  temperature  is,  as  expected,  the  major  source 
of  density  variations,  the  mean  variability  of  salinity  has  been  shown  to  be  responsible  for 
a  significant  equatorward  component  along  the  coast  of  California  and  a  strong  offshore 
component  adjacent  to  Baja.  Based  on  this  analysis,  Batteen  et  al.  (1995)  concluded  that, 
descriptively  and  dynamically,  both  temperature  and  salinity  are  essential  to  accurately 
characterize  the  large-scale  structure  of  the  CCS. 

Here  we  ask  the  question  what  the  combined  wind  and  thermohaline  variability 
effects  are  on  the  large-scale  CCS  circulation.  To  address  this  question  three  model 
experiments,  one  with  climatological  thermohaline  forcing,  one  with  climatological  wind 
forcing,  and  one  with  both  climatological  wind  and  thermohaline  forcing,  are  run  and 
compared  with  each  other.  This  study  is  organized  as  follows:  The  model  and  the 
experimental  conditions  used  in  the  basic  study  are  presented  in  section  2,  while  results  of 
the  three  model  experiments  are  described  in  section  3. 
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n.  MODEL  DESCRIPTION 


A.  MODEL  EQUATIONS 

The  numerical  model  in  this  study  was  originally  used  for  a  coarse  resolution, 
closed  basin  by  Haney  (1974),  and  later  adapted  for  eddy-resolving,  limited  EEC  regions 
with  open  borders  on  the  northern,  western,  and  southern  boundaries  by  Batteen  et  al. 
(1989,  1996)  and  Batteen  (1997).  The  limited  area  EBC  model  is  multi-level,  uses  non- 
adiabatic  primitive  equations  on  a  beta-plane,  and  has  both  baroclinic  and  barotropic 
velocity  components.  The  model  is  based  on  the  hydrostatic,  Boussinesq,  and  rigid  lid 
approximations.  The  governing  equations  are  as  follows: 


du  dp  4  d^u 

(1) 

dt  p^dy  "  A' 

(2) 

o 

II 

(3) 

1 

II 

(4) 

(5) 

(6) 

dS  „4  d^S 

dt  "  “  dz^ 

(7) 

In  the  above  equations,  t  is  time,  (x,  y,  z)  is  a  right-handed  Cartesian  coordinate 

system  with  x  pointing  toward  shore,  y  alongshore,  and  z  upward.  The  corresponding 
velocity  components  are  (m,  v,  w),  T  is  temperature,  S  is  salinity,  p  is  density,  and  p  is 
pressure.  Table  1  provides  a  list  of  other  symbols  found  in  the  model  equations,  as  well  as 
values  of  constants  used  throughout  the  study. 
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For  the  finite  differencing,  a  space-staggered  B-scheme  (Arakawa  and  Lamb, 
1977)  is  used  in  the  horizontal.  Batteen  and  Han  (1981)  have  shown  that  this  scheme  is 
appropriate  when  the  grid  spacing  is  approximately  on  the  same  order  as,  or  less  than,  the 
Rossby  radius  of  deformation,  which  meets  the  criteria  of  this  study.  The  horizontal  grid 
spacing  is  1 1  km  in  the  alongshore  direction  and  8  km  in  the  cross-shore  direction,  while 
the  internal  Rossby  radius  of  deformation  is  ~  30  km.  In  the  vertical,  the  10  layers  are 
separated  by  constant  z-levels  of  13, 46, 98,  182,  316,  529,  870,  1416,  2283,  and  3656  m. 
This  spacing  scheme  concentrates  more  on  the  upper,  dynamically  active  part  of  the 
ocean,  above  the  thermocline. 

The  model  domain  (Figure  2)  is  a  rectangular  region  encompassing  the  west  coast 
of  the  United  States,  from  ~  35°N  to  47.5°N  (1408  km  alongshore),  and  from  ~120°W  to 
132.5°W  (1024  km  cross-shore).  The  coastal  boundaries  of  the  model  domain  are  closed, 
and  have  both  the  tangential  and  normal  components  of  velocity  set  to  zero.  Bottom 
topography  has  been  omitted  to  focus  on  the  roles  played  by  wind  and  thermohaline 
forcing.  The  constant  depth  used  in  the  model  is  4500  m.  A  modified  version  of  the 
radiation  boundary  conditions  of  Camerlengo  and  O'Brien  (1980)  is  used  for  the  open 
ocean  domain  boundaries  to  the  north,  south,  and  west.  Some  spatial  smoothing  is  applied 
in  the  vicinity  of  the  open  boundaries. 

The  model  uses  biharmonic  lateral  heat  and  momentum  diffusion  with  the  same 
choice  of  coefficients  (i.e.,  2.0  x  lO^^cm'^s**)  as  in  Batteen  et  al.  (1989).  Holland  (1978) 
showed  that  the  highly  scale-selective  biharmonic  diffusion  acts  predominantly  on 
submesoscales,  while  Holland  and  Batteen  (1986)  found  that  baroclinic  mesoscale 
processes  can  be  damped  by  Laplacian  lateral  heat  diffusion.  As  a  result,  the  use  of 
biharmonic  lateral  diffusion  should  allow  mesoscale  eddy  generation  via  barotropic 
(horizontal  shear)  and/or  baroclinic  (vertical  shear)  instability  mechanisms.  As  in  Batteen 
et  al.  (1989),  weak  (0.5  cm^s'*)  vertical  eddy  viscosities  and  conductivity  are  used. 
Bottom  stress  is  parameterized  by  a  simplified  quadratic  drag  law  (Weatherly,  1972),  as 
in  Batteen  et  al.  (1989). 
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B.  METHOD  OF  SOLUTION 


Equations  (1)  through  (7)  comprise  a  closed  system  of  seven  scalar  equations  and 
seven  unknowns,  u,  v,  w,  p,  p,  T,  and  S.  The  variables,  u,  v,  T,  and  S  are  prognostic 
variables  whose  time  rates  of  change  are  predicted  from  (1),  (3),  (6)  and  (7),  respectively. 
Although  the  diagnostic  variables  w,  p,  and  p  can  be  determined  from  (3),  (4),  and  (5), 
respectively,  there  are  additional  constraints  imposed  on  p  and  w  by  the  choice  of  the 
rigid  lid  boundary  conditions.  The  vertically  integrated  pressure  can  no  longer  be  obtained 
by  integrating  the  hydrostatic  equation  (4)  for  the  free  surface.  Further,  the  vertically 


integrated  horizontal  velocity  is  constrained  to  be  non-divergent,  i.e.. 


(8) 


which  is  obtained  by  integrating  (3)  and  applying  the  vertical  boundary  conditions. 

For  any  quantity  q,  let  its  vertical  average  be  denoted  by  q  and  its  departure 
(vertical  shear)  by  q' .  From  (8)  the  vertical  mean  flow  can  then  be  described  by  a 
streamfiinction  yr ,  such  that: 


_  1  ^ 

~~ H  dy 


(9) 


I  dy/ 


(10) 


The  streamfiinction  y/  is  predicted  from  the  vorticity  equation,  which  is  derived 


from  the  vertical  average  of  (1)  and  (2).  Applying  the  curl  operator  to  the  vertical  average 
(1)  and  (2),  and  using  (9)  and  (10),  the  vorticity  equation  is 
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d  I  ( (9V  I  1  ( dy/  dH  ^  dy/  dH 


dt  dt^Hydx^  )  H\dy^  )  dx  dx  ^  ^  J 

=jAfx^V— f— — 1 

dx\H  dy  )  dy\H  dx]_  (11) 

where  G  and  F  represent  ±e  collected  contributions  of  the  nonlinear  and  viscous  terms  in 
equations  (1)  and  (2). 

The  vorticity  equation  (11)  is  solved  by  obtaining  an  updated  value  of  ^  by 

application  of  the  leapfrog  (or  every  1 1  time  steps,  the  Euler-backward)  time-differencing 
scheme.  The  associated  value  of  y/  can  then  be  obtained  from: 


1  ( d^y/^  1  ( d^yr^  dyr  dH~'^  dyf  <?//'* 


H\dx^  ]  H\dy 


^  dx  dx  ^  dy  dy 


which  is  an  elliptic  equation.  A  solution  to  (12)  is  fully  prescribed  by  specifying  the 
values  of  y/  on  the  open  and  closed  boundaries  of  the  model  domain.  Currently,  to  solve 

(12),  the  model  uses  an  elliptic  solver  when  there  are  no  variations  in  coastline  geometry 
and/or  topography,  and  successive  over-relaxation  techniques  when  there  are  variations  in 
coastline  geometry  and/or  topography. 

The  vertical  shear  current  («',  v")  is  predicted  from  (1)  and  (2)  after  subtracting  the 


vertical  mean.  The  results  are: 


dt  Po  dx  dz  PoH 


dv'  ~ldp'  A  d^v'  —  T* 

dt  Po  dy  dz  PoH  ^4) 

In  (13)  and  (14),  p',  which  represents  the  departure  of  the  pressure  from  the 
vertical  average,  is ,  using  (4),  expressed  in  terms  of  p  as: 
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/'O 

fO  J  ro 

P'  =j  Pgdt  (J  pgdt)dz 


where  e  is  a  dummy  variable  representing  the  vertical  coordinate. 

The  method  of  solution  consists  of  predicting  V^yr,y/,u',v',T,  and  S  from  (11), 


(12),  (13),  (14),  (6)  and  (7),  respectively.  The  total  current  is  then  obtained  by  adding  the 
vertical  shear  part  to  the  vertical  average  part,  after  the  latter  is  obtained  from  y/  using 


(9)  and  (10).  The  diagnostics  p,  w,  and  p'  are  obtained  from  (5),  (8),  and  (15), 
respectively. 


C.  FORCING  CONDITIONS 


In  this  study,  to  explore  the  effects  of  thermohaline  forcing  on  the  CCS,  seasonal 
temperature  and  salinity  climatological  conditions  from  Levitus  et  al.  (1994)  and  Levitus 
and  Boyer  (1994)  are  used  to  initialize  the  model,  and,  once  a  day,  to  force  the  model  at 
the  western  boundary.  The  seasonal  temperature  and  salinity  forcing  conditions  for  the 
upper  seven  levels,  which  are  initially  assumed  to  be  zonally  homogeneous,  are  shown  in 
Figure  3,  for  the  northern  (47.5  °N)  and  southern  (35  °N)  boundaries  of  the  model 
domain.  Since  the  lower  three  levels  do  not  exhibit  much  horizontal  variation,  they  are 
assumed  to  be  constant  for  each  level.  The  temperature  values  used  for  levels  8  to  10  are 
2.56°C,  2.08°C,  and  2.00°C.  The  salinity  constant  used  for  the  lower  three  levels  is  34.7. 

While  the  temperatures  to  the  south  are  warmer  than  those  to  the  north  at  all  seven 
levels,  only  the  upper  level  temperature  conditions  (Figure  3a  and  3b)  show  significant 
seasonal  variability  with  a  temperature  maximum  in  September  and  a  temperature 
minimum  in  February  throughout  the  whole  region.  Below  these  depths,  (Figures  3c-3g), 
both  the  seasonal  temperature  fluctuations  and  the  temperature  gradient  weaken,  as 
expected.  In  contrast  the  salinity  conditions  at  all  seven  levels  (Figure  3),  which  show  less 
(more)  saline  water  to  the  north  (south),  have  no  significant  seasonal  cycle. 
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To  explore  the  effects  of  wind  forcing,  the  model  is  forced  with  wind  fields  on  a 
2.5°  by  2.5°  grid  from  the  European  Center  for  Medium-Range  Weather  Forecasts 
(ECMWF)  near-surface  wind  analyses  (Trenberth  et  al.,  1990).  The  monthly  mean 
stresses  based  on  twice  daily  wind  analyses  from  1980-1989  have  been  interpolated 
spatially  to  the  8  by  1 1  km  model  resolution  and  temporally  to  daily  wind  values. 

The  wind  forcing  employed  is  displayed  in  Figure  2,  which  depicts  the  seasonal 
winds  starting  with  January  15  (day  15).  The  atmospheric  pressure  pattern  for  January 
(Figure  2a)  has  a  low  (i.e.,  the  Aleutian  Low)  to  the  north  and  a  high  (i.e.,  the  North 
Pacific  High)  to  the  south,  which  results  in  a  wind  divergence  near  40°N.  This  pattern  of 
poleward  winds  north  of  40°N  and  equatorward  winds  to  the  south  continues  through 
February  and  March.  During  April  (Figure  2b)  and  May  the  divergence  in  the  wind  field 
migrates  poleward.  By  June  an  equatorward  component  in  the  wind  field  is  observed 
along  the  entire  domain.  The  strongest  equatorward  winds  are  discernible  from  July 
(Figure  2c)  through  August.  By  October  (Figure  2d)  the  winds  start  to  weaken  throughout 
the  domain,  and  divergence  in  the  wind  field  is  observed  in  the  north.  This  divergent 
wind  pattern  continues  through  November.  By  December  the  wind  divergence  has 
returned  to  ~40°N. 

D.  EXPERIMENTAL  DESIGN 


The  design  of  the  model  experiments  is  as  follows.  Experiment  1  examines  the 
model  response  to  seasonal  temperature  and  salinity  forcing  using  the  seasonal 
thermohaline  forcing  conditions  shown  in  Figure  3.  The  model  integrations  start  from  a 
state  of  rest  and,  once  a  day,  the  model  is  updated  for  temperature  and  salinity  at  the 
western  boundary. 

In  Experiment  2  the  model  is  forced  from  rest  with  seasonal  ECMWF  winds.  The 
initial  mean  stratification  used  is  an  exponential  temperature  profile  with  a  vertical  length 
scale  of  h=450  m.  The  exact  form  is 
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T(z)  =  TB+ATez^ 


(16) 


The  approximation  assumes  Tg  =  2°C  to  be  the  temperature  at  great  depth.  AT  =  13°C  is 
the  increase  in  temperature  between  the  bottom  of  the  ocean  and  the  surface.  This 
temperature  profile  is  the  same  profile  used  by  Batteen  (1989,  1997)  and  Batteen  et  al. 
(1989)  and  was  derived  by  Blumberg  and  Mellor  (1987)  from  available  CCS  observations 
of  the  long-term,  mean  climatological  temperature  stratification  for  the  CCS  region  as  a 
whole. 

To  isolate  the  effects  of  wind  forcing  in  Experiment  2,  the  net  heat  flux  at  the  sea 
surface  is  zero.  In  this  way,  any  heat  flux  that  is  discernible  will  have  been  generated  by 
an  increase  or  decrease  of  sea  surface  temperature  resulting  from  wind-forcing  effects 
(Batteen  et  al.,  1989).  To  accomplish  this,  an  initial  air  temperature  is  chosen  that  forces 
the  net  flux  of  longwave  radiation,  sensible  heat,  and  latent  heat  to  balance  the  heating 
due  to  solar  radiation.  This  air  temperature  is  then  used  in  the  model  for  all  wind  forcing 
experiments.  Any  subsequent  surface  heat  flux  forcing  is  therefore  a  secondary  effect  of 
the  changes  to  the  sea  surface  temperature  due  to  the  wind  forcing. 

Experiment  3  combines  the  thermohaline  forcing  along  the  western  boundary 
(Figure  3)  with  the  seasonal  ECMWF  wind  forcing  to  study  the  model  response  to  both 
types  of  forcing  on  the  CCS.  The  model  integrations  start  from  a  state  of  rest  and,  once  a 
day,  the  model  is  updated  with  ECMWF  winds  and,  at  the  western  boundary,  with 
temperature  and  salinity. 
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in.  RESULTS  FROM  MODEL  SIMULATIONS 


A.  EXPERIMENT  1  -  SEASONAL  TEMPERATURE  AND  SALINITY  FORCING 

1.  Spin-up  Phase 

As  expected,  when  the  temperature  decreases  poleward  (Figure  3),  the  pressure 
gradient  due  to  the  temperature  gradient  establishes  an  onshore  geostrophic  inflow  in  the 
interior  ocean  (Figure  4).  Because  the  initialized  alongshore  temperature  field  is  not 
constant,  the  onshore  flow  varies  between  ~2  and  5  cm/s.  On  approaching  the  eastern 
boundary,  the  onshore  flow  turns  and  forms  a  poleward  boundary  current  (e.g..  Figure 
4a),  which  advects  warm  water  from  the  equatorial  end  of  the  model  domain.  Because  it 
is  continually  augmented  downstream  by  additional  onshore  flow,  the  poleward  current 
increases  in  magnitude  towards  the  pole,  so  that  by  day  45  maximum  velocities  on  the 
order  of  ~30  cm/s  are  observed  near  the  poleward  end  of  the  model  domain  (Figure  4b). 
As  time  progresses,  the  onshore  geostrophic  flow  weakens  in  the  southern  region  to  ~  1 
and  2  cm/s.  The  poleward  boundary  current  also  weakens  so  that  by  day  180  (Figure  4c) 
maximum  speeds  have  decreased  from  ~30  cm/s  at  day  45  to  ~20  cm/s.  A  minimum 
speed  of  ~10  cm/s  is  reached  by  day  255  (Figure  4d),  and  is  maintained  throughout. 

The  onshore  pressure  gradient  due  to  the  thermohaline  gradient  is  sufficient  to 
establish  an  unstable  equatorward  flow  (e.g.,  see  day  63  in  Figure  5).  The  strongest 
equatorward  flow  (~5  cm/s)  is  in  the  coastal,  poleward  end  of  the  domain  and  is 
maintained  throughout. 

Baroclinic/barotropic  instabilities  in  the  poleward  surface  current  and  equatorward 
undercurrent  result  in  the  generation  of  meanders  near  the  coast.  By  day  180  (Figure  4c) 
an  anticyclonic  eddy  develops  at  ~43°N,  subsequently  intensifies  and  propagates 
westward.  From  days  180  to  255  (Figures  4c  and  4d),  other  eddies  develop  equatorward 
of  this  eddy,  and  in  time  intensify  and  propagate  westward. 
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The  surface  (e.g.,  Figure  4c)  and  subsurface  (e.g..  Figure  5)  velocity  fields 
superimposed  on  the  density  fields  illustrate  that  away  from  the  eddy  generation  region, 
the  flow  is  predominantly  parallel  to  the  isopycnals,  as  expected  for  geostrophic  flow.  In 
the  eddy  generation  region,  i.e.,  in  the  poleward  end  of  the  model  domain,  there  is 
significant  advection  of  less  dense  water  offshore  by  the  anticyclonic  eddies. 

2.  Quasi-equilibrium  Phase 

Longer  experimental  runs  show  that  the  system  has  reached  a  quasi-steady  state 
and  that  these  features  continue  to  be  generated  and  maintained.  Using  the  results  of  the 
longer  experimental  mns  (i.e.,  year  3  of  model  simulation  time),  the  model  output  is  time- 
averaged  every  3  days  to  see  the  structure  of  features  in  the  CCS.  As  a  result  of 
geostrophic  inflow  due  to  the  pressure  gradient,  both  a  poleward  surface  current  and  an 
equatorward  undercurrent  are  generated.  Figure  6  shows  the  vertical  stmcture  of  a  coastal 
poleward  surface  current  with  velocities  of  ~8-10  cm/s  overlying  a  coastal  equatorward 
undercurrent  with  velocities  of  -1-6  cm/s  below  ~150-300m  depth  in  the  poleward  part  of 
the  model  domain. 

The  presence  of  both  a  poleward  surface  current  and  an  equatorward  undercurrent 
has  been  observed  by  Hickey  (1979,  1997)  off  the  coast  of  Washington.  During  the  fall, 
the  coastal  poleward  surface  current  has  been  observed  to  have  speeds  of  -5-20  cm/s  and 
to  occur  above  -150-300  m  depth.  The  equatorward  undercurrent  has  also  been  observed 
to  have  speeds  of  less  than  10  cm/s,  and  to  occur  below  -150-300  m  depth.  These 
observations  are  consistent  with  the  results  of  the  fall  model  simulations  near  the  coast  of 
Washington  (e.g..  Figure  6b). 

Due  to  vertical  and  horizontal  shear  instabilities  between  the  poleward  surface 
current  and  the  equatorward  undercurrent,  anticyclonic  eddies  are  generated  in  the  coastal 
poleward  end  of  the  model  domain.  In  time  they  intensify  and  propagate  westward  (e.g.. 
Figure  7).  The  eddies  tend  to  be  0(100)  km  in  size  and  extend  to  depths  of  -300  m  or 
deeper  (not  shown). 
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B.  EXPERIMENT  2  -  SEASONAL  WIND  FORCING 


1.  Spin-up  Phase 

Initially,  the  oceanic  response  to  the  climatological  wind  pattern  (Figure  2)  that  is 
used  to  force  the  model  from  rest,  spins  up  linearly.  In  response  to  the  prevailing 
poleward  winds  in  the  northern  part  of  the  model  domain,  a  surface  poleward  coastal 
current  of  ~3  cm/s  develops  in  the  poleward  end  of  the  model  domain  within  -100  km  of 
the  coast.  After  being  present  for  -45  days,  it  retreats  farther  poleward  following  the 
migration  of  the  Aleutian  Low.  As  expected,  equatorward  wind  forcing  in  the  south 
results  in  an  equatorward  surface  current  (the  CC)  in  the  equatorward  end  of  the  model 
domain.  By  spring,  with  the  arrival  of  the  North  Pacific  High  along  the  west  coast  of 
North  America,  the  current  extends  all  along  the  coast.  A  coastal,  poleward  undercurrent 
(the  CUC)  develops  below  the  surface  equatorward  current,  initially  in  the  equatorward 
end  of  the  model  domain,  and,  during  the  upwelling  season,  along  the  entire  coast. 

As  the  core  of  the  CUC  intensifies,  it  shoals  and  displaces  the  core  of  the  CC 
farther  offshore  (e.g.,  compare  Figures  8a  and  8b).  As  a  result,  there  are  strong  vertical 
and  horizontal  shears  in  the  upper  layers  between  the  CC  and  the  CUC.  Because  this 
current  is  baroclinically  and  barotropically  unstable,  meanders,  filaments,  and  eddies 
subsequently  develop. 

The  equatorward  CC  develops  fairly  uniformly  along  the  coast,  while  the 
upwelling,  observed  as  isolated,  closed  temperature  contours,  is  relatively  patchy  (Figure 
9).  In  particular,  the  coldest,  upwelling  water  is  often  found  at  or  equatorward  of  large 
promontories  (e.g..  Cape  Mendocino  in  Figures  9c  and  9d).  Since  promontories  are  areas 
of  the  irregular  coastline  where  the  alongshore  component  of  the  wind  stress  is  at  a  local 
maximum,  the  CC,  upwelling,  and  growth  of  filaments  should  be  enhanced  in  these 
regions  during  the  upwelling  season  (Batteen,  1997). 
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The  filaments  that  develop  are  “anchored”  off  of  Cape  Blanco,  Cape  Mendocino, 
and  Point  Arena  (e.g..  Figure  9b  and  9c).  The  equatorward  CC  also  form  meanders  in  the 
vicinity  of  capes,  which  intensify  and  develop  into  predominantly  cyclonic  eddies.  In 
time,  these  eddies  coalesce  with  other  cyclonic  eddies  to  form  relatively  large  (-100-300 
km  diameter)  eddies.  After  the  upwelling  season,  the  equatorward  CC  takes  the  form  of  a 
meandering  jet  embedded  with  several  cyclonic  eddies  (e.g.,  Figure  9d). 

2.  Quasi-equilibrium  Phase 

Using  year  3  of  the  results  of  the  longer  experimental  runs,  we  time  average  the 
model  output  every  3  days  for  the  months  of  January,  April,  July,  and  October  to  see  the 
seasonal  structure  of  features  in  the  CCS.  The  results  show  that  in  spring  (Figure  10a)  and 
summer  (Figure  10b),  and  throughout  the  upwelling  season  (Figure  11a),  there  is  a 
coastal,  equatorward  flow  with  speeds  of  -30-80  cm/s  that  leaves  the  coast  between 
-42°N  and  43  °N,  in  the  vicinity  of  Cape  Blanco,  and  meanders  downstream.  The  current 
then  takes  the  form  of  a  meandering  jet.  The  jet  has  alongshore  wavelengths  of  several 
hundred  kilometers  and  cross-shore  excursions  of  several  hundred  kilometers,  and  it  can 
extend  to  depths  of  -700  m  (e.g..  Figure  11b).  South  of  -43°N,  inshore  of  the  jet,  there 
are  both  cyclonic  and  anticyclonic  eddies.  The  cyclonic  eddies  tend  to  form  in  the  vicinity 
of  the  capes.  A  comparison  of  the  spring  (Figure  10a)  and  summer  (Figure  10b)  velocity 
and  temperature  fields  shows  that  as  the  eddies  move  westward,  they  become  embedded 
in  the  meandering  jet.  The  results  of  the  lower  layer  velocity  fields,  averaged  over  the 
upwelling  season  (Figure  1  Ic)  show  that  there  is  a  coastal  poleward  undercurrent,  with 
speeds  of  -10-20  cm/s,  within  -100  km  of  the  coast. 

After  the  upwelling  season,  i.e.,  in  fall  (e.g..  Figure  10c),  a  poleward  flow,  with 
speeds  of  -10-20  cm/s,  develops  within  -100  km  of  the  coast  and  subsequently  replaces 
the  coastal,  equatorward  flow.  Offshore  of  the  poleward  flow  are  several  cyclonic  eddies. 
Farther  offshore,  there  is  a  relatively  strong  (-50  cm/s)  equatorward,  meandering  jet 
embedded  with  both  cyclonic  and  anticyclonic  eddies. 
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In  the  winter  (e.g.,  Figure  lOd),  in  the  poleward  region  of  the  model  domain,  the 
meandering  jet  has  been  displaced  farther  offshore  by  the  westward  propagation  of 
cyclonic  eddies.  In  the  equatorward  part  of  the  model  domain,  the  jet  meanders  closer  to 
shore,  i.e.,  from  ~128°W  at  41°N  to  ~125°W  at  38°N.  Offshore  and  inshore  of  the  jet, 
both  cyclonic  and  anticyclonic  eddies  fill  much  of  the  model  domain.  A  relatively  strong 
surface  poleward  flow  near  the  coast  is  still  present  north  of  ~41°N. 

C.  EXPERIMENT  3  -  SEASONAL  WIND  AND  THERMOHALINE  FORCING 

1.  Spin-up  Phase 

Due  to  the  combination  of  thermohaline  and  wind  forcing,  different  oceanic 
responses  are  expected  depending  on  the  season.  In  the  winter,  in  the  poleward  end  of  the 
model  domain,  the  large  high-to-low  pressure  gradient  due  to  the  warm-to-cold 
temperature  gradient  establishes  an  onshore  geostrophic  flow,  while  the  poleward  wind 
stress  results  in  onshore  Ekman  flow.  On  approaching  the  eastern  boundary,  the  onshore 
flow  turns  and  forms  a  poleward  boundary  current  (e.g..  Figure  12a).  In  the  equatorial  end 
of  the  model  domain,  the  smaller  pressure  gradient  and  the  equatorward  wind  stress 
results  in  weak  onshore  geostrophic  flow,  offshore  Ekman  flow,  and  a  coastal 
equatorward  surface  current  (e.g..  Figure  12a). 

During  the  upwelling  season  (-April  to  September)  (e.g..  Figures  12b-12d),  the 
combination  of  a  weakened  pressure  gradient  and  increased  equatorward  winds  over  the 
entire  model  domain  lead  to  a  strengthening  of  equatorward  flow  all  along  the  coast  and  a 
weakening  (strengthening)  of  onshore  (offshore)  flow.  As  in  Experiment  2,  the  upwelling 
is  relatively  patchy  with  the  coldest  water  found  near  coastal  promontories  (e.g..  Cape 
Blanco  in  Figure  12d).  The  subsurface  structure  of  the  currents  shows  that  there  is  a 
poleward  undercurrent  (e.g..  Figure  13a).  The  undercurrent  is  within  -150  km  of  the 
coast,  extends  from  -150  to  greater  than  700  m  depth,  and  has  a  core  velocity  of  -5  cm/s. 
Typical  core  velocities  for  the  undercurrent  range  from  -5  to  20  cm/s.  The  surface  coastal 
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equatorward  flow  extends  to  ~  150  m  depth-near  shore  to  ~500  m  depth  offshore.  The 
core  of  the  surface  current  is  within  -100  km  of  the  coast  and  has  a  typical  core  velocity 
of  -50-70  cm/s.  The  upwelling  season  model  simulations  are  consistent  with  observed 
depths  and  velocities  of  both  the  surface  equatorward  current  and  the  poleward 
undercurrent  (e.g.,  Huyer  et  ah,  1991). 

In  the  fall,  as  expected,  the  surface  poleward  flow  strengthens  in  the  coastal 
poleward  end  of  the  model  domain  in  response  to  both  the  strengthening  of  the  pressure 
gradient  and  the  return  of  poleward  wind  stress  in  the  region.  The  subsurface  structure  of 
the  currents  shows  that  in  the  poleward  end  of  the  model  domain,  there  is  an  equatorward 
undercurrent  (e.g..  Figure  13b).  The  undercurrent  is  within  -30  km  of  the  coast,  extends 
from  -150-700  m  depth,  and  has  a  core  velocity  of  -2  cm/s.  Typical  velocities  for  the 
undercurrent  range  from  -2  to  5  cm/s.  As  in  Experiment  1,  the  fall  model  simulation  in 
the  poleward  end  of  the  model  domain  is  consistent  with  the  observed  depths  and 
velocities  of  both  the  surface  poleward  current  and  the  equatorward  undercurrent  near  the 
coast  of  Washington  (Hickey,  1979, 1997). 

2.  Quasi<equUibrium  phase 

Longer  run  times  (-3  years)  of  the  model  simulation  show  that  during  the 
upwelling  season  many  of  the  features  simulated  in  Experiment  2,  such  as  meanders, 
eddies,  a  temperature  front,  and  filaments  (e.g..  Figures  10  and  11)  are  also,  as  expected 
with  the  dominance  of  wind  forcing,  simulated  in  Experiment  3.  A  comparison  of  these 
features  for  Experiments  2  and  3  during  the  upwelling  season  (Figures  11a  and  14, 
respectively)  shows  that  the  results  are  qualitatively  similar  except  that  there  is  a  tighter 
temperature  gradient  inshore  of  ~125°W  between  Cape  Mendocino  and  Point  Arena  in 
Experiment  3.  This  must  be  due  to  the  onshore  geostrophic  flow  resulting  from  the 
alongshore  pressure  gradient.  An  examination  of  recent  CCS  observations  by  Strab  and 
James  (1995)  shows  that  during  the  upwelling  season  (e.g.,  see  Figure  2a  of  Strub  and 
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James,  1995)  the  temperature  front  between  Cape  Mendocino  and  Point  Arena  is  found 
inshore  of  ~125°W,  consistent  with  the  model  simulation  of  Experiment  3. 

Due  to  the  thermal  wind  relationship,  stronger  surface  and  subsurface  currents 
inshore  of  ~125°W  in  Experiment  3  are  also  found  between  Cape  Mendocino  and  Point 
Arena.  For  example,  Figures  15a  and  15b  show  that  the  surface  coastal  jet  with  velocities 
of  ~10  cm/s  north  of  Cape  Mendocino  (Figure  15a)  intensifies  to  ~50  cm/s  at  Cape 
Mendocino  (Figure  15b),  while  the  undercurrent  with  velocities  of  ~15  cm/s  north  of 
Cape  Mendocino  intensifies  to  ~25  cm/s  at  Cape  Mendocino. 

In  the  poleward  end  of  the  model  domain,  a  comparison  of  subsurface  currents 
(not  shown)  for  Experiments  2  and  3  shows  that  the  velocities  for  the  imdercurrent 
decrease  from  ~45  cm/s  in  Experiment  2  to  ~20  cm/s  in  Experiment  3.  The  reduction  in 
the  velocities  for  the  poleward  undercurrent  is  likely  due  to  the  net  result  of  a  poleward 
undercurrent  generated  by  wind  forcing  opposing  an  equatorward  undercurrent  generated 
by  thermohaline  forcing  in  this  region.  These  lower  undercurrent  velocities  are  more 
consistent  with  available  observations  of  the  undercurrent  (e.g.,  Huyer  et  a/,.  1991). 

Horizontal  maps  of  the  upper  layer  mean  kinetic  energy  (MKE)  and  eddy  kinetic 
energy  (EKE),  averaged  over  the  duration  of  the  upwelling  season,  are  shown  in  Figures 
16a  and  16b.  Maps  of  MKE  and  EKE  are  suggestive  of  where  the  mean  and  eddy  energy 
sources  are  to  be  found  (Holland  et  al,  1983).  A  comparison  of  Figures  14  and  16  shows 
that  high  values  of  MKE  and  EKE  are  found  all  along  the  coastal  and  offshore  axes  of  the 
equatorward  jet,  and  in  the  offshore  regions  south  of  -42°N. 

A  comparison  of  Figures  16a  and  16b  shows  that  maximum  values  of  MKE  and 
EKE  occur  in  the  same  region,  i.e.,  inshore  of  ~125°W  between  Cape  Mendocino  and 
Point  Arena.  Note  that  the  MKE  values  are  larger  than  the  EKE  values  in  the  region.  This 
is  consistent  with  the  results  of  Batteen  (1997),  which  showed  that  the  eddies  are 
generated  from  instabilities  of  the  mean  equatorward  current  and  the  poleward 
undercurrent  via  baroclinic  and/or  barotropic  instability  processes. 
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Maps  of  EKE  for  each  month  during  the  upwelling  season  (not  shown)  show  that 
the  maximum  EKE  always  occurs  between  ~37  to  39°N  and  between  124  and  128°W. 
The  maps  also  show  that  high  values  of  EKE  are  present  throughout  the  upwelling  season 
between  Cape  Mendocino  and  Point  Arena,  while  off  of  Cape  Blanco,  EKE  values 
steadily  increase  from  small  to  large  values  as  the  upwelling  season  progresses,  with 
highest  values  occurring  in  July  and  August.  These  results  suggest  that  the  greatest  eddy 
activity  off  of  Cape  Blanco  occurs  after  June,  while  eddy  activity  between  Cape 
Mendocino  and  Point  Arena  occurs  regularly  throughout  the  upwelling  season. 

Since  Kelly  et  al.  (1997)  have  most  recently  used  observational  results  to 
investigate  the  variability  of  near  surface  EKE  in  the  CCS,  we  qualitatively  compare  the 
EKE  results  of  our  model  simulation  with  their  observations.  Both  Kelly  et  al.  (1997)  and 
our  results  show  the  location  of  the  maximum  EKE  to  be  between  ~36  and  40°N,  and  at 
~125°W  during  the  upwelling  season.  This  region  of  maximum  EKE  coincides  with 
increased  equatorward  flow  in  both  our  model  simulations  and  their  observations.  Both 
Kelly  et  al.  (1997)  and  our  results  (not  shown)  also  show  that,  south  of  ~40°N,  on  a 
seasonal  time  scale  the  region  of  monthly  maximum  EKE  migrates  westward  to  ~128°W. 

After  the  upwelling  season,  i.e.,  in  the  fall,  the  combination  of  wind  and 
thermohaline  forcing  results  in  the  seasonal  changes  in  the  currents  in  the  poleward  end 
of  the  model  domain.  In  particular,  the  reversal  of  winds  from  equatorward  to  poleward 
combined  with  the  intensification  of  the  alongshore  pressure  gradient  results  in  a  shoaling 
of  the  poleward  undercurrent  near  the  coast  of  Oregon  and  Washington,  as  seen  in  the 
latitudinal  sequence  of  Figures  17a- 17c.  These  results  are  consistent  with  the  hypothesis 
of  Hickey  (1989,  1997)  that  a  surface  poleward  current  can  result  from  a  shoaling  of  the 
undercurrent. 

On  the  basis  of  these  results  we  can  conclude  that,  even  though  the  effects  of  wind 
forcing  dominate  the  CCS,  the  additional  effects  of  the  thermohaline  forcing  result  in  the 
following:  the  development  of  a  poleward  surface  current  and  an  equatorward 
undercurrent  in  the  poleward  end  of  the  model  region;  an  onshore  geostrophic 
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component,  which  results  in  a  temperature  front  and  stronger  surface  and  subsurface 
currents  between  Cape  Mendocino  and  Point  Arena;  and  a  region  of  maximum  EKE 
inshore  of  ~125°W  between  Cape  Mendocino  and  Point  Arena,  associated  with  the 
temperature  front.  These  results  seem  credible  given  that  the  phenomenological  behavior 
of  the  model  simulation  has  been  shown  to  be  qualitatively  similar  to  recent  large-scale 
hydrographic,  altimetric,  drifter  and  moored  observations  of  the  CCS. 
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Longitude  (degrees  west) 


Figure  1.  Generalized  circulation  schematic  of  the  classical  climatological  California 
Current  System  (CCS).  The  broad,  slow  surface  equatorward  California  Current 
(CC)  overlies  the  poleward  California  Undercurrent  (CUC)  along  with  the  Inshore 
Countercurrent  (IC),  known  as  the  Davidson  Current  (DC)  north  of  Point 
Conception  and  as  the  Southern  California  Eddy  (SCE)  south  of  Point  Conception. 


20 


Day  15 


Longitude  (degrees) 

20  m/s - > 


Figure  2.  Domain  and  climatological  winds  of  the  model  for  the  California  Current 
System  (CCS)  off  the  Western  United  States.  The  domain  of  the  model  is  bounded 
by  35°N  to  47.5°N,  120°W  to  132.5°W.  The  cUmatological  (1980-1989)  ECMWF 
winds  used  in  Experiments  2  and  3  in  m/s  are  shown  here  for  (a)  January,  (b)  April, 
(c)  July,  and  (d)  October.  Maximum  wind  vector  is  20  m/s. 
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13  Meter  Thermohaline  Forcing  At  132.5W 


Figure  3.  Time  series  plot  of  monthly  temperature  and  salinity  fields  used  as 
seasonal  forcing  in  Experiments  1  and  3.  The  **’  symbol  represents  data  at  35°N, 
132.5W,  while  the  ‘+’  plot  represents  data  at  47.5°N,  132.5W  for  levels:  (a)  13  m,  (b) 
46  m,  (c)  98  m,  (d)  182  m,  (e)  316  m,  (f)  529  m,  and  (g)  870  m  depth. 
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46  Meter  Thermohaline  Forcing  At  132.5W 


20  h 


O 

Vl5[ 

3 
03 

U. 

(1) 

0) 
t- 


35N 


* 


3K  M 

+  47.5N 
+  +  + 


SK 


+  + 


6  8 
January-December 


10 


12 


34.5 


_  34 

3 

M 

CL 

533.5 

'c 

To 
CO 


33  h 


32.5  L 
0 


*35N 

*  *  * 
+  47.5N 


*  * 


xc  * 


6  8 
January-December 


10 


12 


b 


26 


Salinity(psu) 


98  Meter  Thermohaline  Forcing  At  1 32. 5W 


January-December 


January-December 

c 


27 


182  Meter  Thermohaline  Forcing  At  132.5W 
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316  Meter  Thermohaline  Forcing  At  132.5W 


i 

*35N 

- - 

- - 

1 

+  47.5N 

* 

*  3K 

3K 

* 

51^ 

* 

- 

+ 

+  + 

_ 1 _ 

+ 

+ 

+ 

+ 

-1  — 

+ 

+ 

+ 

+ 

_ _ 

- 

20 

15 

10 


6  8 
January-December 


10 


12 


34.5 


32.5 


6  8 
January-December 


e 


29 


529  Meter  Thermohaline  Forcing  At  132.5W 
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DEPTH  :  46m 
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u  .  V - >  50.0  LONGITUDE 

Figure  4.  Density  contours  and  velocity  vectors  at  46  m  depth  for  Experiment  1  at 
days  (a)  3,  (b)  45,  (c)  180,  and  (d)  255.  In  all  the  velocity  fields  presented,  to  avoid 
clutter,  velocity  vectors  are  plotted  every  third  (fourth)  grid  point  in  the  cross-shore 

(alongshore)  direction.  Contour  interval  is  0.2  gm/cm^;  maYimnm  velocity  vector  is 
50  cm/s. 
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DEPTH  ;  529m 
T  (DAY)  :  63 
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Figure  5.  Density  contour  and  velocity  vectors  at  529  m  depth  for  Experiment  1,  day 
63.  Contour  interval  is  0.2  gm/cm^J  maximum  velocity  vector  is  10  cm/s. 
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Figure  6.  Cross-shore  section  of  meridional  velocity  (d)  at  46°N  for  Experiment  1 
during  the  third  year  of  model  simulation,  time-averaged  for  days  (a)  1  -  363  and  (b) 
270-363.  Contour  interval  is  1  cm/s. 
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Figure  9.  Temperature  contours  and  velocity  vectors  at  46  m  depth  for  Experiment 
2  at  days  (a)  225,  (b)  255,  (c)  270,  and  (d)  315.  Here  and  the  following  figures,  CB 
refers  to  Cape  Blanco,  CM  to  Cape  Mendocino,  and  PA  to  Point  Arena.  Contour 
interval  is  1°C;  maximum  velocity  vector  is  100  cm/s. 
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Figure  11.  Time-averaged  plots  for  the  upwelling  season  for  Experiment  2  of  (a) 
temperature  contours  and  velocity  vectors  at  46  m  depth,  (b)  cross-section  of  i)  at 
4U°N,  and  (c)  velocity  vectors  at  316  m  depth.  Contour  interval  is  1°C  in  Figure 
11a  and  5  cm/s  in  Figure  11b.  Maximum  velocity  vector  is  100  cm/s  in  Figure  11a 
and  80  cm/s  in  Figure  11c. 
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Figure  12.  Density  contours  and  velocity  vectors  at  46  m  depth  for  Experiment  3  at 

days  (a)  45,  (b)  87,  (c)  180,  and  (d)  255.  Contour  interval  is  0.2  gm/cm^;  maYimuin 
velocity  vector  is  50  cm/s. 
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Figure  13.  Cross-shore  of  v  in  the  coastal  region  for  Experiment  3,  days  (a)  180  at 
39°N  and  (b)  300  at  46°N.  To  avoid  clutter,  only  velocities  up  to  10  cm/s  are  plotted 
for  Figure  13b.  Contour  interval  is  5  cm/s  for  Figurel3a  and  1  cm/s  for  Figure  13b. 
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Figure  15.  Cross-shore  sections  for  Experiment  3  in  the  third  year  of  model 
simulation,  time-averaged  for  days  90-270  at  (a)  41'’N  and  (b)  40.2°N.  Contour 
interval  is  5°C. 
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Figure  16.  Horizontal  maps  at  13  m  depth  of  (a)  mean  kinetic  energy  (MKE),  and 
(b)  eddy  kinetic  energy  (EKE)  for  Experiment  3  in  the  third  year  of  model 
simulation,  time-averaged  for  days  90-270. 
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Figure  17.  Cross-shore  sections  for  Experiment  3  in  the  third  year  of  model 
simulation,  time-averaged  for  days  270-363  at  (a)  37°N,  (b)  43°N,  and  (c)  46°N. 
Contour  interval  is  5  cm/s. 
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Table  1.  Values  of  Constants  Used  in  the  Model 


Constant 

Value 

Deflnition 

To 

278.2°K 

Constant  Reference  Temperature 

^0 

34J 

Constant  Reference  Salinity 

Po 

Density  of  Sea  Water  At  Tq  and  Sq 

a 

BBE^sam 

Thermal  Expansion  Coefficient 

P 

7.5  X  10-^ 

Saline  Expansion  Coefficient 

K 

10 

Number  of  Levels  In  Vertical 

Ax 

Cross-Shore  Grid  Spacing 

Ay 

Alongshore  Grid  Spacing 

H 

ESBBSIBSHHI 

Total  Ocean  Depth 

800  s 

Time  Step 

fo 

Mean  Coriolis  Parameter 

8 

Acceleration  of  Gravity 

Am 

2  X  10^’ cm^ 

Biharmonic  Momentum  Diffusion  Coefficient 

Ah 

2  X  10^^  cm'^  s'* 

Biharmonic  Heat  Diffusion  Coefficient 

Km 

0.5  cm^  s'* 

Vertical  Eddy  Viscosity 

Kh 

0.5  cm^  s'* 

Vertical  Eddy  Conductivity 
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